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ABSTRACT 

This paper describes the development of mathematical models 
which reflect the effects of tracking system measurement errors (i.e., 
noise and bias) and tracking system location uncertainties on RMS er- 
rors in spacecraft position and velocity. Through the application of the 
method of weighted least squares, a generalized mathematical model is 
developed which permits simulation of tracking the spacecraft with 
several tracking stations in the simultaneous or nonsimultaneous track- 
ing mode. In Part I only simultaneous tracking was considered. 

Mathematical models are also developed for: (1) effects of a priori 
knowledge on the RMS errors in spacecraft position and velocity; (2) the 
ellipsoid of errors in spacecraft position and velocity; and (3) target 
body impact parameters, the errors in impact position and impact time. 

The last two models are especially important because they pro- 
vide the techniques to use when describing the geometry for the RMS 
errors in spacecraft position and velocity at discrete points along the 
trajectory and at any targeting or impact point. Finally, a numerical 
example is given which serves to illustrate the effects of known error 
sources on the RMS errors in spacecraft position and velocity. 
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TRACKING SYSTEMS, THEIR MATHEMATICAL 
MODELS AND THEIR ERRORS 
PART II: LEAST SOUARES TREATMENT 


by 

W. D. Kahn and F. O. Vonbun 
Goddard Space Flight Center 


INTRODUCTION 

The rapid increase in the number of analyses of space missions during the past few years 
and in the increasing complexity of the missions has intensified the need for more accurate de- 
termination of the position and velocity of spacecraft. A realistic estimate of the errors in any 
such determination should include errors resulting from a variety of sources. These are the er- 
rors associated with the basic measurements, the uncertainty of the station locations, the basic 
physical constants, and the mathematical models used to describe the orbit. In addition, random 
errors, which are errors subject to statistical analysis, must be distinguished from bias errors, 
which are errors that stay constant or change slowly during the measurement process. 

This report presents a detailed error analysis utilizing a weighted least squares approach, 
and shows how all of the types of errors mentioned above can be combined for different tracking 
systems. This is important, from the practical point of view, since on each space mission many 
tracking systems and stations are being utilized to determine the spacecraft trajectory. In order 
to design and evaluate a tracking system properly, all significant error sources mentioned must 
be considered. 

Most of the formulas are taken from Part I of this paper (Reference 1) keeping changes at a 
minimum. Many of the equations given in Part I are included for the reader T s convenience. Par- 
ticular attention is paid to the bias errors in the measurements and in the locations of the stations. 
This paper shows that the bias errors are predominant when only short tracking times can be 
utilized for the trajectory determination. To cite a practical example, when a spacecraft is in- 
serted into a parking orbit, as in the case of the Apollo Missions, only a tracking ship located in 
the Atlantic Ocean can "see" the insertion and then only for 2 to 4 minutes (Reference 2). This 
time is too short to evaluate bias errors and to correct for them properly. Evaluation of some of 
these bias errors is possible when more than one station is utilized for the trajectory determination 
and when time is less important. This is the case for multiple earth orbits, transfer orbits to the 
moon and planets, and orbits with long periods (days rather than hours). 

Errors, for proper interpretation, are treated in a geometric sense; that is, the error ellip- 
soids for position and velocity are calculated. The principal axes of these ellipsoids are determined 
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and their projection along a space trajectory is treated in three dimensional form. The hyperel- 
lipsoid treatment (6 -dimensional) has been avoided since it lacks geometric interpretation, 
whereas the error ellipsoids of the position and velocity can easily be visualized and transformed 
into impact miss parameters at the moon or planets. 


ANALYSIS 

Part I of this paper was primarily concerned with the development of the error equations for 
the fundamental spacecraft tracking systems. These systems are: (1) Radar Systems which meas- 
ure range, azimuth, and elevation; (2) Interferometer Systems which measure direction cosine and 
direction cosine rates; (3) Ranging Systems which measure range only; and (4) Range Rate Systems 
which measure Doppler only. These are called fundamental systems because they can be combined 
to handle all known present and future tracking systems. As an example, consider the USBS,* 
which is presently used in the Apollo program. This system measures range, range rate, 
and x-y angles, which are readily converted to azimuth and elevation (Reference 3). This 
system can be presented as a combination of a radar system and a range rate system. 
Previously the determination of RMS errors in spacecraft position and velocity was made 
by assuming k different tracking stations making one measurement each at the discrete 
points along the trajectory which were visible to each of the k tracking stations. In brief, simul- 
taneous observations were necessary to determine the spacecraft position and velocity errors 
without utilizing the orbit itself (Reference 4). This method is important when extremely accurate 
spacecraft position and velocity are to be obtained. SECOR* is a typical example of a system 
which does not make use of the orbit and thus is not subject to the errors introduced in the equa- 
tions of motion. 

For the analysis presented here, error sources are of secondary importance and major em- 
phasis is given to the errors themselves which influence the trajectory determination. Although 
the errors and their sources are interrelated, a separation is assumed in order to reduce the size 
of Part II. The position and velocity errors are determined by combining the measurements, n k , 
made by each of the k tracking stations during the time the spacecraft is "radio-visible" from the 
station. "Radio-visibility" means that the station position vector intersects the spacecraft antenna 
at sufficiently close range to provide a signal strength within given limits; e.g., 3 db, 10 db, etc. 
(Reference 5). The solution for these errors is obtained by the weighted least squares method, 
combining all of the measurements for the stations viewing the spacecraft either simultaneously 
or nonsimultaneously. The basic requirement here is that the number of measurements has to 
exceed the number of unknown parameters to be determined. In most cases these unknown 
parameters are the six orbital parameters or the position and velocity vectors. Since only the 
errors in the state vector (6x1 vector composed of the spacecraft position and velocity components) 


*USBS stands for Unified S-Band System which is a system combining range and range rate measurements with voice and telemetry and 
television. 

tSECOR is the Sequential Electronic Correlation Ranging System. 
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are solved for, a spacecraft need not be visible to each of the k tracking stations at the same time. 
Obviously there are many more observations than unknowns to be solved for. 

Random as well as bias errors are treated in this paper. In its colloquial use, the word "ran- 
dom" is applied to any method of choice which lacks aim or purpose. In the true sense, a measure- 
ment subject to random errors must be repeated many times in order to ascertain the most prob- 
able value and to assure that any future measurement lies within a specified interval of confidence. 
This is to some extent in conflict with what is being done in spacecraft tracking, since no quantity 
considered can be measured more than once because the spacecraft is in motion. In most practi- 
cal cases, however, the relative changes in the magnitude of these quantities are small during 
the measurement interval so that proven techniques for error analysis can be applied. 

Bias or systematic errors are in general those which stay constant from measurement to 
measurement or fluctuate in periods which are larger than necessary to make the measurements 
to determine a quantity within given limits. Throughout Part Et, bias errors are not solved for and 
thus add to the spacecraft total position and velocity errors. However, bias errors in the measure- 
ments and tracking station locations are considered in the overall analysis along with the tracking 
system’s uncertainties. For single station solutions, bias errors play a larger role than does 
measurement noise. As an example, if only measurement noise is used in the analysis, the errors 
in the state vector tend to approach zero; this situation is obviously meaningless in the practical 
sense. Incorporation of the bias error effects, which cannot be avoided, will in practice wipe out 
this phenomenon automatically and add considerably to the errors in the state vector (Reference 6). 
The error equations in this paper are derived from the variational equations given in Part I. 

Since Part n treats all random measurements in a least squares sense, the letters a, /3, 7, and 
s (representing the type of tracking systems in Part I) are no longer applicable in the same way. 
The reason is that in a practical sense each system makes a different number of measurements. 

To be as consistent as possible, these letters are now used as subscripts indicating the type of 
tracking systems used in the analysis; that is 

a designates radar, 
ft designates interferometer, 
y designates ranging system, and 
S designates range rate system. 

The total number of measurements made by each of these systems are designated as follows: 


A. Definition of Symbols 

The following symbols will be used throughout this paper. 
4> - Geodetic latitude of the tracking station. 


3 


L 



k - Geodetic longitude of the tracking station, 
h - Height above geoid of the tracking station. 

a 0 - Equatorial radius of spheroid used to represent the Earth, (For 
= 6378.388 km). 

e* - Square of eccentricity for spheroid, (for Hayford Spheroid = 

e G - Greenwich sidereal time at 0 h Universal Time (U.T.), obtained 
almanacs. 

<9, - Greenwich sidereal time at U.T. of observation, 
t - U.T. of observation, 
a - Azimuth of the object being tracked, 
e - Elevation of the object being tracked, 
r - Slant range of the object being tracked. 
a s - Semi-major axis of satellite orbit. 
k' - Right ascension of satellite. 

S - Declination of satellite. 
e s - Eccentricity of satellite orbit, 
i - Inclination of satellite orbit. 

s 

- Longitude of ascending node. 
oj s - Argument of perigee. 

E s - Eccentric anomaly. 

G s - True anomaly. 

P - Period of revolution of satellite. 

s 

N - Radius of curvature along prime vertical. 
p s - Magnitude of the radius vector to the satellite in the inertial coordinate system. 
h s - Satellite altitude above the Earth or other central body, 
j v | - Magnitude of velocity vector of the satellite. 
y s - Flight path angle. 

B. Error Analysis for Spacecraft Position and Velocity 

Mathematical Models 

The variational equations corresponding to the fundamental tracking systems are presented 
in this section. These equations differ from those derived in Part I in that they reflect the meas- 
urement variation with respect to the bias in the measurement. For a single measurement in the 


Hayford Spheroid a 0 

.0067226700223). 
from the various 
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inertial coordinate system these equations are as follows:* 


1. Radar Systems (Measuring r , a, e) 

D (3i a x3) SK (3X1) =L (3i a x3) 5 X ( 3 x 1 ) “ (j TR )(3i a x3) S S ( 3 x 1 ) + D (3i a x3) S ^(3xl) 


where 


= Bias in Measurement 


and 


i = 1, 2, . . . , n . 


2. Angle and Angular Rate System (Interferometer— Measuring a, e, a, e) 

5K (3xl) =L ° (2i „ 3) SX (3X1) ~ (F a J TR \2 i(3 x 3 ) SS (3X1) +D “ (2i x3) Sp 


and 


D° s i 3 = M° 5X (3 n + L° 8X - (F V J t r) 5S.J j. + D° Sp 

a (2i J gX3) ( 3X1 ) a (2i fl x 3 ) C3X1) a (2i fi x 3 ) ( 3xl > Xa 3 > ^ ) a (2i^<3) ( 3xl ) 


where 


l^g - 1, 2, . . . , ^ 


3. Ranging Systems (Measuring r only) 


%„> SK <-> -% y „> - (F . J, “'<v> ss <-> 


where 


i — 1 » 2, 3, •••> n^ 


4. Range Rate Systems (Measuring f only) 

D r 8K (3 xn =«! 8X 3 + L° s X 3xn -(F r VJ T R) x3 SS +DJ Sp 

r (igX3) ( 3X 1) r ( i gX 3 ) (3X1) r ( i gX 3 ) ( 3xl > r ( 1 S X 3 > ( 3X1 ) r (igX 3 ) 


•All pertinent matrices used in the equations discussed in this paper are fully defined in Appendix a. 


L 


(B-l) 


(B-2) 


(B-3) 


(B-4) 


(B-5) 
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where 


= 1, 2, ■ 


Error Equations in Combined Form 

Present day tracking systems can be presented as combinations of the so-called basic sys- 
tems which are used in this paper. The mathematical models of the error equations for such sys- 
tems can therefore be formed by writing in partitioned matrix form combinations of Equations 
(B-l) through (B-5). A solution for the RMS error in the state vector (position and velocity) can 
then be made by applying the method of weighted least squares. Listed below are well known com- 
binations of advanced types of present day tracking systems. In all instances, they represent 
single measurements. 

1. Range and Range Rate System (Measuring r , r, a, e, or equivalent) 

This system (for example, the USBS and the Goddard Range and Range Rate System) is a com- 
bination of the models shown in Equations (B-l) through (B-5). The variational equation in parti- 
tioned matrix form is: 


) 


(DMO (3x1 


L (3x3) 

0 (3x3) 

M° 

L° 

r (lx3, 

_ f ( 1 x 3 ) 




‘{3xl) 


(4x 1) 


(4x 6 ) (6-1) 


(-J t R) 


(3x3) 


<- F r*J T R) (lx3 ) 


ss, 


(3xl) 


D (3«3) 


0 (1*3) 

D? 

(1x3, 


fe P(3*l 


(4x3) 


_^C3*n_ 

(4x6) (6xi) 


(B-6) 


In essence, Equation (B-6) is applicable for other combinations of systems, such as (1) a radar 
system and a Doppler system and (2) systems measuring range only or range rate only with an 
interferometer system measuring a and 6 or their equivalent. 


2. Range Rate - Interferometer System (Measuring f , a, e, or equivalent) 

This combination of systems can be mathematically represented by combining Equation (B-2) 
and (B-5) as follows: 




L° 

8 ( 2 x 3 ) 

°(2x3) 


fX (3»l) 


-(f 8 J T B\ 2 x 3 ) 

>S (.«.) + 

D° 

a ( 2 x 3 ) 

0 ( 2 x 3 ) 


bP (3xl) 



M° 

r (l*3) 

L° 

r (lx 3 )_ 


i x <3..,j 


_-(F r VjTR) (lx3L 

°(lx 3 ) 

D° 

(1*3) 


_^P(3x!)j 
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3. The Angle and Angular Rate System or Interferometer (Measuring a, e , a, and k or equivalent) 
Matrix Equations (B-2) and (B-3) are combined: 


<nS»» ( 2.« 


L° 

*(2* 3) 

°(2x 3 ) 


" 8X C3*.> 


-ff.J T R\ 2 x3) 

S S(3x„ + 

D° 

“(2x3) 

0 (2x 3 ) 



UP 

L (2*3) 

*(2x3) 


_ SX (3«1)_ 



®'2x3) 

D° 

"( 2 xar 

(4xl) 


(4x6) (6* 1 ) 

(4x3) 


(4 


§P (3*1) | 
S P(3xi) 


(B-8) 


Formulation of the Variance- Covariance Matrix for the Spacecraft 1 s State Vector 

The fundamental problem considered here is that of estimating the RMS errors in the state 
vector when uncertainties in the tracking system measurements, measurement biases, and bias 
errors in station locations are known. It should be pointed out that station location errors (except 
for ships and aircraft) are to be considered as constant (or bias errors) for orbit determination, 
since the station location does not change as the spacecraft transits a tracking station. For il- 
lustrative purposes, the mathematical model for the range and range rate system is used to de- 
termine the errors in the state vector resulting from utilizing such a tracking system. In order 
to facilitate future matrix manipulations, Equation (B-6) is rewritten in abbreviated form as 
follows: 


( 4 xl) 


A (4x 6 ) X (6xl) + B (4 x 3) ®(3xl) + ^(4x 6 ) ^(6x1)** 


(B-9) 


where 


(4x1) 


(D §K) ( 3 x j ) 
s »(i«D 


(4xl) 


(4x6) 


(6xl) 


i 

(3X3) ; 

°(3x 3 ) 

M° 1 

_ ' ( 2 x 3 ) i 

L° 

r (lx 3 )_ 

' SX (3 X1) ' 


«(3M,_ 

(6xl) 


(4x6) 


♦Matrix A is not to be confused with the generalized Matrix A given in Part I. 
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B 


(4x3) 


J T ^(3x3) 
[C-FVJ t R) ( 1x3) _ 


U (4x3) 


^(3xl) C (3xi) 


8s . 


Ss , 


( 3x 1 ) 


v 4 x 6 ) 


®(3x3) ; 

°(1*3) 

0 (1*3) j 

D° 

r (l*3)_ 


— I (4x6) 


and 


^(6xl ) 


"^>(3X0 


_(SP)( 3X i) J 


(6xl) 


Equation (B-9) is a linearized equation (terms of a Taylor expansion terminated with the first 
order partial derivatives) which states that the deviations in the measurements represented by the 
column matrix Y (4xl) are dependent on the deviations in spacecraft position and velocity represented 
by the column matrix X (6xl) , the deviations in the location of the instrument making the measure- 
ment represented by the column matrix S (3xl) , and the bias errors represented by the column 
matrix p . A solution of the matrix Equation (B-9) will now be made utilizing the method of 
weighted least squares. Before doing so, several preliminary conditions must be met; these are: 

1. Proper evaluations in spacecraft position and velocity over the interval (t 0 , t ). 

2. Proper weighting of the measurements in accordance with known accuracies associated with 
these measurements. This means giving greatest weight to the most accurate measurement and 
the least weight to the least accurate measurement. 

To fulfill condition 1, it is necessary to find a set of variables which are independent of time. 

Since, in Equation (B-9), the variation in the state vector X (6xl) is time dependent, this 
equation in its present form cannot be solved in the least squares sense. In order to overcome 
this, X (6xl) has to be expressed as a linear combination of another set of variables. In general, 
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this can be represented by the following equation: 


(6x1) - 


^X(t) 

3X(t 0 ) 


%> = t 0 )] 

W f G. V 1 \ U 


(6 X 1) 


(6x6) U (6xl) 


J (6x6) 


(B-10)* 


where t ) //r _ is the "state transition matrix” which linearly transforms the state of the 

0(6x6) 

vector X (6xl) . This transformation is valid provided that linearity is preserved in the interval 

(to.*). 

The elements of the state transition matrix are the first order partial derivatives of the state 
vector components at time t with respect to time t 0 . For a precision trajectory these partial 
derivatives are determined by numerical integration of the equations of motion. 

A precision trajectory is often approximated with a patched conic trajectory (Reference sec- 
tion entitled "Target Body Impact Parameters"). Using a patched conic trajectory for error 
analysis permits expressing the state transition matrix elements in closed form. This approach 
significantly reduces computation time on large digital computers and usually provides a good 
estimate to the actual state transition matrix. 

In this paper the elements of the state transition matrix are expressed in closed form using 
the Keplerian elements as the bridge between the state at time t and at time t 0 . The state transi- 
tion matrix using Keplerian elements is derived as follows: 


(B-ll) 


*(6xl) - ^^^(6x6) A (6xl) 


X °(6xl) ~ P(t 0 ) (6x 6 ) A ( 6 x 1 ) ^ 


where 

X /£ 1N = Variation in state evaluated at time t, 

(6x1) 

X n = Variation in state evaluated at time , 

(6x1) 0 


P(t) 


(6x6) " 


{■ 


3(x,, x 2 , x 3 , 
B(Q , i , co , 



(6x6) 


*Dr. S. F. Schmidt of Philco Corporation - Western Development Laboratories pointed out to the authors the need to introduce the state 
transition matrix into Equation (B~9) in order that the unknown variables comprising the vector X (6xl) can be evaluated in the least squares 
sense. 
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and 


A ( 6 * 1 ) = 


sn 

s 

Si 

s 

$00 

s 

Sa 

s 

Se 

s 

- S 7 s- 


, Variation in the Keple rian Elements . 


(6x1) 


As can be seen, the matrix P(t) (6x6) has to be the Jacobian matrix transforming the variation 
in the orbital elements into the variation in the state vector. The matrix P(t), is fully derived 

(6xg) 

in Appendix a of this paper. From Equations (B-ll), the following result is obtained: 


'(6 


<1) - [P(t) P 1 (^^(6x6) X 0 (6xl) ~ ( Pf V)( 6x6 


w (6xl) 


[<&(t, 


t 0^ (6> 


6) X 0 


(6x1) 


(B-12) 


The matrix product [PP^ 1 ] (6x6) represents the state transition matrix* in the interval (t 0 , t). 
It is now possible to solve for the errors in the state vector in the least squares sense by in- 
troducing (B-12) into (B-9), resulting in: 


V (4xl) ^4*6) X 0 (6 x 1) + 


D (5 

(4x3) °(3xl) 


+ ^(4x4) P(4x 1 ) 


(B-13) 


where 


*"(4x6) ^PP Q )^ 4x6) 

The unbiased linear estimate of the unknown vector x o is obtained from an abbreviated form 
of Equation (B-13): (6 * 1} 


(4xl) 


r X 

' (4x 6 ) *0 


(6xl) 


(B-14) 


Equation (B-14) represents the linear terms of a Taylor T s expansion about the deviations which 
comprise the column matrix X (higher order terms are neglected, see Part I, Appendix b). 

76x11 

This equation has to be modified in order to reflect the omission of higher order terms of the 
Taylor’s expansion as well as the effect of the bias errors and station location uncertainties. The 


*W. H.Goodyear of IBM Federal Systems Division has derived the state transition matrix in terms of Cartesian coordinates. This form 
of the transition matrix bypasses numerical difficulties which arise whenever the spacecraft’s trajectory is non-elliptic; i.e., para- 
bolic or hyperbolic. For a description and derivation of the state transition in Cartesian form see Reference 7. 
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modified expression for y (4xl) then is: 

Y (4xl) = F (4x6) X 0 (6xl) + C (4xl) ' (B-15) 

Matrix T <4xl) is a column matrix of the residuals resulting from the linearization process. 

The matrix Equation (B-15) represents a single observation made by the range and range rate 
system. In that equation, six components of the variation in the state vector are estimated in the 
least squares sense. It is therefore necessary for the number of measurements, n, to be much 
larger than the number of unknowns. In terms of Equation (B-15), a "single measurement" means 
a measurement of one range, one azimuth, one elevation, and one range rate. Therefore, n meas- 
urements of such a system actually means 4n single observations. Therefore, Equation (B-15) is 
rewritten: 


1 ( 4n x 1 ) 


' ( 4n x 6 ) ^0 


(6x1) 


( 4n x 1 ) 


(B-16) 


where 


4n > 6 . 

Equation (B-16) represents the general equation of a range and range rate tracking system 
making n measurements (this is the same as a radar and a Doppler). All errors in the system are 
attributed to random noise. Extending Equation (B-16) to two tracking systems necessitates the 
following modification: 



x 1) 


r, 

*(4nj«6) 

V „ + 

C 1 

( 4n j x 1 ) 

Lv, 

Xl) 


r 2 

(4n 2 x6) 

(6x1) 

£ 2 

(4n 2 M)_ 


4(nj+n 2 )xl 4(n 1 +n 2 )x6 4( ni +n 2 )xl 


(B-17) 


The general variational equation for i tracking stations, i > 2, is obtained by expanding the system 
shown for two stations by Equation (B-17) as follows: 


Y. 


V^njxe) 


£ 1 

1 (4n 1 x i ) 

y 2 

(4n 2 «l) 

= 

r 2 

(4n 2 x6) 

X + 

(6x1) 

'V/ 

£ 2 

(4n 2 xl) 

Y. • 

*(4n.«l)_ 


r. 

*(4n i x 6 ) 


1 (4n | «X)_ 


(B-18) 


11 



Or, abbreviated form. 


\j 'V# 

K ( 4 -£xl) = G (4-tx6) X 0 (6xl) + £ (4-txl) 


(B-19) 


where 


is Y i - 1 

k = l 

and 


't >>> 6 


The Gaussian condition for solving a system of equations is that: 

0= (* T W _1 *) (1x1) =Min* (B-20) 

= (K T -%G T ) (lx4/t) W^ 4x44) (K-GX 0 ) (Ux1) = Min , (B-21) 


where is the covariance matrix associated with the measurements. This matrix weights 

measurements in accordance with their relative accuracy, so that the greatest weight is given to 
those measurements having the smallest uncertainties associated with them. 

The condition that </>= minimum can be obtained by differentiating Equation (B-21) with re- 

spect to Xn , that is: 

^(6*1) 


AJ 


= - 2 (g t r 1 K) (#x1) + 2 (G t r 1 G) (#x6) \ 0(6xi) = 0 


(6X1) 


(6x1) 


(B-22) 


The minimum condition also requires that [3 2 $/3Xg] (6x6) be positive definite. 
From (B-22), 

aV 

3X 2 


(G T W" 1 G) 


(6x6) 


♦See Reference 1, Equation (B-6) Page 29* 
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If (G T W -1 G) (6x6) is positive definite, (eigenvalues are all positive) then </> is a minimum. From 
(B-22) the estimate to the column matrix X 0 is: 


\ 6xti = (G T W' 1 ^(6X6) W_1 *>(6.1) • (B_23) 

Let a matrix N be defined to represent only the random noise in the measurements. Then if 


and 


E(KK t ) (44xU) =E(NN t ) (4XxU) =w (4 ^ x4<t) 


(B-24)* 


measurements from different tracking stations are assumed to be statistically independent (un- 
correlated). The following expansion of matrix W (4 ^ x4 ^ ) reflects this assumption by having null 
matrices off the main diagonal. 


Vf 


(4-tx4-t) 


( 4 n j* 4 n j ) 

o * 

0 


(4n 2 x 4n 2 ) 

0 


0 

0 

3 (4n * 4 n ) 


w. 


1 ( 4 n.x 4 n i ) 


( 4 / t x 4 ^) 


The covariance matrix in terms of spacecraft position and velocity is derived by using 
Equations (B-23) and (B-24) as follows: 


E «0*S\6-6 ) =( C(rW " 1 G)7l 


6 ) G (6 


< 4 ^) ^( 4 ^x 4 ^) E ^^ T ^( 4 ^x 4 ^) 4 ^) G ( 4 "£x 6 ) 


(G T W' 


1 GV 1 
w; (6> 


6 ) 


*E( ): The expectation operator (Reference 8). 
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= (G T *r‘ G)-J x6) G[ 6xU) E(NN T ) (Ux4/e) W^ xU) G^.,, (G^W' 1 G)"*, 


6 ) 


(G W G) ( 6 * 6 ) G*.^ W(4^ x4 ^ ) ^(4 '£x4'C) ^(4 / £x4'£) ^(4^x6) 


(G T W' 


G) 


“1 

(6x6) 


(g t w _ 


lG >(6x6)^ T W- 1 G) (6x6) (G T r 


G \ 6 x 6 ) 


(G T W 


'i Q\-l 
° ; (6> 


6 ) 


(B-25) 


In terms of i tracking stations, Equation 25 when written in expanded form using equation (B-18) 
is: 


e ^ 0 XJ) (6x6) = 


-i-i 


L E kk w 


k=1 j k =1 


(6x6) 


(B-26) 


where 


n k > 6 - 


k > 1 , 
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and 



To obtain the biased estimate for the range and range rate system, let the error in the meas- 
urement reflect the effect of measurement noise, measurement bias, and station location uncer- 
tainties as follows: 


+ B (4^x3i)®(3i xl) +^(4^x41) "(4ixl) 


(B-27) 


where 
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By substituting Equation (B-27) in (B-23), the biased estimate to the variations in the state vector is: 


'v _ ^ ^ 4 'V 'y r \y 'V/ 

\ 6 «d =(G w G) <«- 6 ) {g w < N + BS + Cb ^>(« 


*D 


(B-28) 


The covariance matrix in the state vector for the biased estimate then is:* 


E <X 0 Xj) (6x6) = (G t W -1 G) ( j x6) G[ 6xU) *-!t x4<) (E(NN t ) + BE(SS t ) B t + CE(b b T ) C T } (UxU) g ( 4 ^ x6) (G T w 

Let it be assumed that: 


^*Vf 'V* rp <\y <\p X 

E(NS ) = E(Nb ) = E(Sb T ) =0 


then 


E <%>X T o\e« 6) = <G r *~ 1 V~ ( l,> <W) ^U.4-e, tw + BE(SS T ) B T +CE(bb T ) C t } (44xU) W-U x 4 ^G (Ux6) (G^G)^, 


6 ) 


(G T r i B) (6X 3 i) E(ss T ) (3 . x31) (B T r i G) ( 3. x6i) (G T r i G^ x6) 


+ <G T r G );‘ x6) (G T r j c) (6x4i) E(bb T ) (4ix4i) (is T w Ji 6> <4i . e) dfir 1 ©-;.,, 


(B-29) 


or, in expanded form, 


e(\ 0 x;) = 


i 

E 


(6x6) k =i 


K 





i 

E 




iw = 1 





i 

E 

k = l 


k 

E <r ‘*‘‘ B ‘ > -. 


iv = 1 




k 'i k 


ik-i 


E 

k*l 






r >\ 


U* 1 


E 

k=l 


T <? V E>; 


4—1 

ifc"! 


Jk 



E E 


j k =1 




K 

2] ^ Tw " r ^i, 

ik =1 



E 

k = l 


k 

E 

i k =1 


<c*:‘ V,. 


V*1 


(B-30) 


‘ *This proof was suggested by W. Goodyear, IBM Federal Systems Division. 



where 


and 


[ E ( S S\] ( 3«3) = < 


(3x3) 


k (4x4) 


r 

V 2 o 0 0 


< 

O 

prL 

o 

o 

> 


0 0^0 



0 0 0 v 2 ' 


V. 

T 

k j( 


-I K J(4x4) 


where 


k = index for the k th tracking station. 


The errors in the state vector as obtained from Equation (B-29) or (B-30) are those which 
correspond to a least squares solution in which the effects of the bias errors were included but not 
solved for. This has the effect of obtaining a nonoptimal least squares estimate of the errors 
in the state vector. To obtain the optimum least squares estimate of the errors in the state vector 
requires solving for corrections to the tracking station location as well as for corrections to the 
measurements. The optimal least squares estimator is now to be derived for a single range and 
range rate tracking system. For n measurements Equation (B-13) reads: 

Vn*l) = r (4nx6) ^0 f6xl j + ®(4nx3) ^( 3 x 1 ) + ^(4nx4) ^(4x1) • (B"31) 

Assuming that the deviations in station location and in the measurements remain constant 
necessitates introducing two equations of constraint, which are stated as follows: 


S °(3«l) 


(For Deviations in Station Location) j 
(C Po\4n*i) = c (4nx 4) P(4«i)( For Deviations in the Measurements^ 


(B-32) 
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From a combination of Equation (B-31) with (B-32) the following matrix equation results: 


(4nx 1) 


u (3xl) 


(U.1 


( 8n + 3 ) x 1 


— 

®(4nx3) 

— 

f (4nx6) 

^(4n x4) 

_ _ 

^(3x6) 

I (3«3) 

0 (3x4) 




4n x 6) 

0 (4nx3) 

r 

( 4n x 4 ) 


(8n + 3)x 13 


( 6 * 1 ) 


(3x1) 


^(4x1) 


(13x1) 


(B-33) 


or, in abbreviated form, 


I - P T 

*"*( 8n + 3) x l ~~ cf( 8n+3 ) x 1 3 


(B-34) 


The omission of higher order terms resulting from the linearization process necessitates restating 
Equation (B-34) as follows: 


^ 'X/ 

^(8n + 3)xl = §(8n+3)xl3 ^ (13xl) +£ (4n+3)xl 


(B-35) 


As described in Equations (B-20) through (B-23), a best estimate for the unknown column 
matrix T is given by the following: 




(B-36) 


where 



E (YY T ) ( 4 n x4„) 

0 

0 

W E(L Q L 0 ) C g n+3)x ( 8n + 3 ) - 

0 

E (S 0 ^(3x3) 

0 


0 

0 

{C E(pp cn (4nx4n) _ 


(8n+3)x (8n +3) 


and 


E(VY T ) (4nx4n) =W (4nx4n) 

< E ^ 0 s;) ( 3«3) = W (3*3, 

^ E (P 0 P o\4nx4n) = tW ^(4nx4n) 
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It then follows that the optimal estimate for the errors in the components of the column matrix 

~ . 

T is: 


where 


e <t 0 tJ) (1 3 X1 3, = (§ t w- 1 §);; 3x1 3> 


(B-37) 


1 

E(XoXj) 

E^S 7 ) 


E ^O^o\l3Kl3) 

E(SXj) 

E(SS t ) 

E(S£ T ) 


E(frxJ) 

E(pS T ) 

E(?? T > 


(§ T W _1 §\; 3*13) 


n 

^2 (r T w _1 n. 


22 (B T ?r» n. 

j=i 


22 ( cT w_i 


n 

22 0" T w _1 B)j 


|22 w_i B >j + w_i 


n 

22 (C T W _I B) j 


22 ( rT W_1 c >j 

j=i 


(B t W" 1 C). 


n 

22 (c x w_i Oj + 


1 (13x 13) 


By inverting the matrix on the right hand side of Equation (B-37), the optimal estimate of the co- 
variance matrix for the state vector E(X 0 Xj) is: 


e cw (6 * 6) = 

where 


” [ A I1 ~ A 13 A 33 A 13 “ ^12 “ A 13 A 33 ^3^ ^ 22 A 23 A 33 A 23^ 1 ^2 A 23 A 33 ^ 
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(6*3) 


J 2 (rr wi B) i 

j*i 


and 


A ,3 3 


n 

(r T W' 1 C). 

j = l 


(6x4) 


A - A T 
A 21 “ A 12 


a 22 


(B t W ' 1 B). + w' 1 


1 (3*3) 


a 23 s \tl (BTflrlc) i[ 

Ife J(3X 4) 


A - A T 

a 31 ~ 13 ’ 


A - A t 

32 “ ^23 * 


A 33 ~ 


(C T W 1 C) j + 


' (4x4) 


u (4v 4 ) 


0 0 0 

r cos f 0 0 

0 r 0 

0 0 1 


(4x4) 


The nonoptimai estimate of the errors in the state vector for a single tracking station ob- 
tained by setting index i of Equation (B-30) to one is: 


E ^\ T \ 6 x 6) = 


{ A Ii + A U A !2 » a I 2 M + \s A [s A 7i} 


(6x6) 


(B-39) 
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A numerical evaluation of Equations (B-38) and (B-39) will show the estimate for the errors in 
the state vector to be low when Equation (B-38) is used and high when Equation (B-39) is used. The 
most reasonable estimate of the errors in the state vector probably is obtained from a combination 
of the non-optimal estimator with the optimal estimator. This can be accomplished by solving for 
some but not all of the bias errors in the least squares sense. 

Effects of a Priori Knowledge on the Errors in the State Vector 

As the powered flight portion of a trajectory is terminated and the free flight portion is in- 
itiated, some knowledge about the RMS errors in spacecraft position and velocity at the transition 
point can be assumed. Such knowledge is defined here to be a priori knowledge for the free-flight 
trajectory; that is to say, the following RMS errors are assumed to be known at the time the space- 
craft is injected into orbit. 

r} h - RMS error in spacecraft height above the earth, 

7] v - RMS error in speed, 

s 

rj - RMS error in flight path angle, 

T} a - RMS error in injection azimuth, 

- RMS error in right ascension (equivalent to RMS error in longitude), and 
^ - RMS error in declination (equivalent to RMS error in latitude). 

In order to include the effects of these RMS errors in the RMS error in spacecraft position 
and velocity, the following procedure is used. The transformation relating [sh , Sv s , 8y s , Sa., 8 8(8)] 
with [ox j, Sx 2 , 8 x 3 , 8x 1? 8x 2 , 8x 3 ] is given by the following matrix equation: 


X(t) = y (t) p(t) , (B-40) 

(6x1) (6x6) (6x1) 


where 


( 6 * 1 ) 


sx 


(3*1) 


sx, 


(3*1) 


(6*1) 


*(0 ^ 
(6*6) 


~J 

3(x r x 2 , x 3 , ij, i 2 , i 3 )] 

f 

1 

L 3<H lWs) a it V, 8) J 

t 


21 



and 


P(6*l) 


'Sh 
S v 

s 

S ^ s 

8a. 

i 

BK 9 

8 ( 8 ) 


(6xl) 


The time that a spacecraft is inserted into orbit corresponds to the epoch time t. Equation (B-4G) 
evaluated at time r reads: 


XW(6«1) = ^ (T) (6x6) P W (6,1) 


(B-41) 


By using the result from Equation (B-12), Equation (B-41) can be restated as follows. 


XW(6.|, = t P ^^ 1] (6,6)\ 6xl) = ^(6*1) 

P^(6xl) = ^t ( 6s1) = [P T P 0 1 ] (6x6) *0 (6xl) 


(B-42) 


or 


B = Jy" 1 [ 0 ( t , t )] X \ 0 

Pt (6x1) \ T (6x6) 0 (6*6)J V* 1 ) 


Equation (B-42) is combined with Equation (B-19) as follows: 


'V 

K 


1 

VO 

K 

V 

rr 

(S* 

1 

'X- 

v 

€ (4^x 1) 

}r_ 


_{▼"> *(T, t 0 )} (6xlL 

X Q + 

°(6xl) 

_°(6xl)_ 


(4^+6)*l (4^+6)*6 (4^+6)xl 


(B-43) 


The covariance matrix reflecting the effect of a priori knowledge of the state vector com- 
ponents is derived in accordance with Equation (B-25) as 




(B-44) 
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where 


V> = [ p MP _1 (t 0 )] (6 , 6) • 

For the nonoptimal least squares estimator, the influence of a priori knowledge introduces a modi- 
fication of Equation (B-39) as follows: 


e (W ( 6 «6> = [B Ii + B n \2*\ T 2 B Ii + B il *13 b;;] 


where 


l (6x6 


} s [v +» T {t t E(p r p r T ) Y^r 1 «] ^ 


6x6) 


(B-45) 


Similarly, using matrix [B xl ] (6x6) in place of matrix [A u ] (6x6) in Equation (B-38) will incorporate 
the influence of a priori knowledge into the optimal least squares estimator. To incorporate the 
effects of a priori knowledge into the covariance matrix of spacecraft position and velocity, it is 
sometimes presumed that the RMS errors are uncorrelated; that is, 


E (I| T >(6*6) = TraCe -P T >(6 x6) = 


Tjjj 0 0 0 0 o" 

0 7)1 0 0 0 0 

s 

0 0 tj* 0 0 0 

/ s 

0 0 0 r? o 0 

i 

0 0 0 0 1)2, o 

0 0 0 0 0 Tjf 


(B-46) 


By neglecting the covariance elements in the covariance matrix, a conservative estimate for 
the error volume is always obtained (References 9 and 10). Whenever the covariances in the 
variables (h, v s , y s , a., 5) are available, they are to be included in the covariance matrix. 

This means that E( pp T ) <6x6) is nondiagonal and symmetric; or, in other words: 




n, % 


(B-47) 
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where 


'ij 


= Correlation Coefficient 


f- 1 < p. . < 1 for i / j 
(j°ij =1 f°r i = j 


?7- _ = Covariance. 


The effect of a priori knowledge is very significant whenever the amount of tracking data is 
limited. As the amount of tracking data is increased, the influence of the a priori knowledge di- 
minishes to a point of insignificance. 


The Ellipsoid of Errors in Spacecraft Position and Velocity 
The normal distribution function for the vector (Reference 11) 


(6*1) 


6X 


(3xl) 


6X 


(3x1) 


(6x1) 


8x, 


Sx„ 


Sx ( 


Sx r 


L_ -3J 


(6xl) 


is given by 


f(8x lT Sx 2 , . . .8x 3 ) 


(2 n) 3 j A | 1/2 


exp - ^ (X T A” 1 X) (1; 


l) ’ 


(B-48) 


where 


A ( 6x6) “ E(XX T )( 6 x 6 ) (The variance-covariance matrix in 
spacecraft position and velocity) 


and 


| a| s Determinant of matrix a. 
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r 


The normal distribution function (B-48) is constant for points of the hyper ellipsoid (X T &‘ 1 X') (lxl) = k 
where k is a constant; that is to say, the probability level that a random point (Sx 1? Sx 2 , • • ■ Sx 3 ) 
will fall inside the hyperellipsoid is dictated by the value assigned to the constant k. It is more 
logical to consider the error ellipsoids in position and velocity separately, primarily because 
greater physical significance can be associated with such surfaces rather than with the hyper- 
ellipsoid. In so doing let 


SX(3 Xl)S (X-X) (3xl > 

^( 3 X 3 ) s E(SXSX T ) (3x3) (Covariance matrix in spacecraft position), 

SX (3xl) ,(X-*) (3xl) . 


in which 


and 


A ( 3 x 3) ~ E (SX $X t ) (3x3) (Covariance matrix in spacecraft velocity), 


X = Mean value of X 


X = Mean value of X , 


|a| = Determinant of matrix & 


I A [ = Determinant of matrix A 


The normal distribution functions corresponding to position and velocity vectorial components 
are as follows: 


f(Sx r Sx 2 , 8 x 3 ) 


3/2 | T I 1/2 


(2ttY /£ I A 


exp - 1/2 (SX T A ” 1 SX) (lxl) 


25 



and 


f(Sx lT Sx 2 , Sx 3 ) 


(277)3/2 I A 1 1/2 


exp 


- [f (si^SX)^. (B-49) 


Each of the normal distribution functions (B-49) is constant for points on the ellipsoids: 


(SX t a _1 sx> (Ixl) = k, 


(B-50) 


and 


( 5 X T A ' 1 SX) (lxl) = fc, , 


where and k 1 are constants. Therefore, the probability that a random point (Sx 17 Sx 2 , Sx 3 ) will 
fall inside the error ellipsoid for position is (Reference 11): 


p = 



t 2 e“ l /2 dt , 


(B-51) 


so that when k 1 = 1.5382, p = .5. The corresponding ellipsoid is called the 50 percent error 
ellipsoid. The same probability level is associated with the velocity error ellipsoid. 

The probability levels for different values of k x , as generated by Equation (B-51), are given by 
the following table (Reference 11): 


1.101 

.25 

1.538 

.50 

2.027 

.75 

2.500 

.90 

2.795 

.95 

3.368 

.99 


To simplify the mathematical analysis of the shape and orientation of the error ellipsoids 
(B-50), a technique known as ' ’Transformation to the principal axis” is used. Described math- 
ematically, the equation of an ellipsoid falls into the general category of the ’’quadratic form.” 

It can then be stated that there exist linear transformations 

SX (3xl) = T (3*3) ^(3x1) (B-52) 
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which reduce the quadratic form 


(5X T A" 1 SX) (lxl) =k x 


into a linear combination of squares: 


te T S *]( 1 . 1) = k 1 


(B-53) 


(B-54) 


where 


D s (T T a" 1 T) 

U (3k 3) “ 1 1 3 x 3 ) 


Matrix D (3x3) is diagonal; that is, the elements of the matrix have the following properties: 


d.. 

‘ j 



Vi 

i = j 


Elements d. . are often referred to as the eigenvalues of the quadratic equation and the columns 
of matrix T (3x3) in (B-52) as the eigenvectors. Thus the orientation of the error ellipsoid can be 
established with reference to a designated coordinate system. 

To show the transformation to the principal axis in more detail, let Equation (B-53) be re- 
written as follows: 


3 3 

y ' c. . Sx. Sx. = k 1 , 
/ . u i j i 

j- 1 i=l 



(B-55) 


where 


^ij -*(3x3)y A (3 x3) 


and 


c. . 
1 j 


c. . 

J 1 
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It remains to reduce Equation (B-55) to the form 


3 



i=l 


(B-56) 


where 


[d..] 


ii (3 k 3) 


(T t A ' 1 T) 

(3x 3 ) 


To do so, let it be assumed that element c 12 (= c 21 ) is the largest off-diagonal element of 
matrix A" 1 , which is to be eliminated by rotation of axes through an angle 6 in the x x x 2 -plane. 
The transformation equations to be used are: 


SXj 5= Xj - Xj = ^ COS 0-^2 s *~ n @ 

S X 2 = x 2 - X 2 = s i n 0 + co s 6 > • 

Sx 3 =X3-X3=f 3 

Substituting Equation (B-57) into (B-55), the coefficient of <f 2 is: 

2 cos 6 sin 6 (c u - c 22 ) + 2c 12 (cos 2 0 - sin 2 0 ) . 

Eliminating the term (i.e., by setting the coefficient of € x % 2 equal to zero), 

2C 12 

tan 26 = — — - 

< C 22 - C U> 


(B-57) 


(B-58) 


(B-59) 


then defines the rotation through the angle 6 . 

Thus if element c 12 is the largest off-diagonal element in matrix a' 1 , the linear transforma- 
tion given by Equation (B-52) will reduce Equation (B-55) to the form given by Equation (B-56). 
For this example matrix T is of the form: 


(3x3) 


cos 6 - sin 6 0 

sin 6 cos 6 0 


0 


0 


1 


(3x3) 
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If, on the other hand, Equation (B-59) is: 


2c i3 

tan 29 - — , 

< C 33 " C lP 

then the corresponding Matrix T would be defined in the following manner: 


(3*3) 


cos 9 0 sin 9 

0 1 0 

sin 9 0 cos 9 


J (3x3) 


In general, a symmetric matrix can be diagonalized by a linear transformation such as given 
by (B-52), and the matrix T (3x3) of this transformation equation always will be diagonal. 

In practice the diagonalization process of a symmetric matrix is accomplished by iteration; 
that is, 


(3x3) 


(3x3) 


(3x3) 


a T A' 1 T ) 

1 “ 1 l ; (3x3) 

(T T D T 'i 

' *2 U l 1 2 ; (3x3) > 

(j n t) 

V n n“l n ; (3x3) 


(B-60) 


For each iteration the largest off-diagonal element can be any d. ( i t j ) . Again it is to be men- 
tioned that the elements d.. (i = j ) in matrix D are the eigenvalues and the columns of matrix 
(T = Tj T 2 . . • # T n ) (3x3) are the components of the eigenvectors. 

Target Body Impact Parameters — The Errors in Impact Position and Impact Time 

In most cases, the impact of a target body by a spacecraft is best described by the errors in 
impact position and impact time. A set of parameters developed by W. Kizner (Reference 12), 
which is based on a technique applied in scattering theory in atomic physics, is used to describe 
these impact errors. 

To properly derive the planetary impact parameters, it is first necessary to describe the 
M patched conic trajectory" concept. Although a precision trajectory for a spacecraft is computed 
by either direct integration of the equations of motion (such as Cowell 1 s method) or by a varia- 
tional scheme (such as Encke's method), these precision trajectories can be approximated by a 
patched conic trajectory. The term "conic" arises from the fact that the solution to the two-body 
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problem in celestial mechanics is a conic section, such as a circle, ellipse, parabola, or hyper- 
bola. The type of conic section is dependent on the energy in the system relative to the central 
body. A conic trajectory is always referenced to a particular central force field. When a conic 
approximation to an actual trajectory (initially referenced to the central force field of body 1) 
reaches the region of dominance of the central force field of body 2, it is referenced to body 2’s 
central force field. The effect of switching the field of reference is that the characteristic of the 
conic section is altered. In the case of a patched conic trajectory from the earth to the moon, the 
trajectory usually is an ellipse with respect to the earth’s central force field. Upon entering the 
region of dominance of the moon’s central force field, the trajectory is referenced to the moon. A 
conic section characterized by a hyperbola results. The process of combining the two conics in 
order to approximate an integrated trajectory is known as "patching conics." 

The concept of the patched conic trajectory makes it possible to develop a set of parameters 
which describes the "miss" geometry at a target body such as the moon or planet. Figure 1 de- 
picts the geometry as a spacecraft approaches a target body. The target plane, given by the unit 
vectors R° and f°, contains the impact vector B. This plane is perpendicular to the asymptote of 
the approach hyperbola, which is the conic resulting from switching from the earth’s central force 
field to the force field of the target body. A unit vector s ° is therefore defined as that unit vector 


/ 

/ 



Figure 1— Target body impact parameters. 
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directed along the incoming asymptote which is normal to the impact plane. The vector T° represents 
an arbitrary vector normal to S° (usually assumed parallel to some fundamental plane such as the 
ecliptic, celestial equator, etc.). The vector R° is normal to both s° and t°; i.e., r° = (S°x?°). 

The miss vector B is that vector which is perpendicular from the center of the target body to the 
asymptote of the approach hyperbola. 

Figure 2 depicts in greater detail the IT impact plane” and the miss vector B. In actuality the 
vector B is not the real "miss” vector but that vector which is directed from the center of the 
target planet to the point of penetration® (Figures 1 and 2) of the trajectory hyperbola and the 
"impact plane.” The real impact point is the trajectory penetration point (p*) (Figures 1, 2, and 3). 
From the "miss” vector B projection onto the unit vectors T° and R°, the direction of the "miss” 
with respect to the penetration point® is determined. A more accurate solution for describing 
the direction of "miss” is with respect to the actual penetration point (?). 




The necessary mathematical derivation for 

these parameters is based on the description of Figure 3-Target body-approach hyperbola, 

the "miss” parameters. Therefore, let the 

spacecrafts position and velocity vectors be defined within the target body f s region of dominance. 
The coordinate reference frame in which these two vectors are referenced is chosen for convenience 
(such as target body — equatorial, target body — ecliptic, etc.). The conic trajectory with reference 
to the target body is a hyperbola. In order to derive the "miss” parameters using Figure 3, it is 
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necessary to define the following vectors: 


P° = — — (unit position vector of spacecraft) , 


v° = — — (unit velocity vector of spacecraft) , 


-3 x ^ (unit momentum vector directed 
s s along trajectory rotational axis 

|p x v | and perpendicular to trajectory 
s 5 plane). 


A unit vector directed to the point of closest approach from the target body is defined as: 


where 


p° = cos 6 P° - sin 6 (n° x p° ) 


I - 1 a s ^ e s ~ ^ (polar form for equation 
Ps ~~ 1 + e cos e of a hyperbola) , 


{| vj } 2 = Pj 


2 1 \ (energy integral for 


p a 

~ s 1 s 


hyperbolic orbits) , 


Q = true anomaly or polar angle , 


semi-major axis, 


e = eccentricity , 


p T B = gravitational constant of target body. 


A unit vector perpendicular to po is: 


q° = sin 6 po + cos 0 (n° x pO . (B-6 

S g s s \ 

The unit vectors p° andqo are directed along the positive and 77 axes of Figure 3. Conse- 
quently the unit vector in the direction of the incoming asymptote is: 


cos a p° + sin a qo 


(B-64) 
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where 


(a? + b *)»/* e s ’ 


b „ /, 1 V 

3 1 n a ~ — = II \ 

o l + b s^ 1/2 \ e s/ 


i\l/2 


and 


b = a (e 2 - l) ly/2 

s s v s 

By definition, the "miss” vector § is normal to S°: 

B = | B | sin ap° - | B | cos aqo t (B-65) 


where 


1 B I = {a +a (e - 1)} sin a 

11 s s s 


(e 2 - 1) 


1/2 


The miss vector with respect to the actual penetration point 


is: 


B* = 


sin a p< 


- B* 


cos a q° 


(B-66) 


where 


1 + e s sin a 

Let the arbitrary unit vector T°, which lies parallel to a chosen fundamental plane (i.e., 
equatorial plane, ecliptic plane, etc.) be defined as perpendicular to S°. A unit vector R° then is 
defined as normal to T° and s°, i.e., 


(S° x T°) 


(B-67) 
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The definition of 5 ° implies that the vector R° is normal to the fundamental plane; thus its com- 
ponents are {0, 0, 1} . Therefore it is possible to write the components of vector T° as follows: 


where 


*J» o _ 


S 2 S, 

1 + J 


(si + s p 1 " <s i +s 2 2 y ' 2 


i = { 1 , 0 , 0}, T = {0, 1 , 0}. 


(B-68) 


Using Equation (B-65) or (B-66) in combination with Equations (B-67) and (B-68), it is possible 
to resolve the "miss" vector B or B* into its components and thus obtain a measure of how much 
in-plane and out-of-plane the "miss" of the target point was (referenced to some chosen funda- 
mental plane). That is, 


(B • t°) = in-plane miss distance 


and 


(B * R°) = out-of-plane miss distance . 

Transformation of the covariance matrix in spacecraft position and velocity (at the time of 
impact with a target body) into a covariance matrix in terms of the "miss” parameters furnishes 
more significant information than can be derived from the RMS errors in the state vector 1 s com- 
ponents at impact. To make this transformation, a third variable is defined. This variable can 
either be energy, usually defined as C 3 , or impact time defined at t 0 , where 


c 


3 



(B-69) 


The following transformations are used to express the covariance matrix in terms of the "miss" 
or "impact" parameters: 


n?L s v.> 




B(B • R°) 


> l ax, 


SX„ 


»o / C 6 * 1 ) 

u '(1*6) 
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(B-70) 


and 


SX„ 


( 6 * 1 ) 




( 6 * 1 ) 


If impact time is chosen as the third variable, the covariance matrix in terms of the "miss" 
parameters is expressed as follows:* 


where 


E(M* M* T ) 


(3x3) “ 



e ^oXJ) (6 , 6) 



(B-71) 


M 


* 

(3xi) 


S(B • f°) 
S(B • R°) 


and 


3(5 • T°) 



3(B • R°) 

3X 0 


3t o 

BX 0 


J (3x6) 


♦The partial derivatives in matrix (BM* x are in practice evaluated numerically and not in closed analytic form due to their 

complexity. Numerical evaluation of these partial derivatives is also more efficient on high speed digital computers. 
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C. Discussion and Interpretation of Results 


The mathematical theory described in the previous sections was used in developing the "Erran 
Error Propagation Computer Program (Reference 13). This computer program was principally 
used for studying the effects of errors inherent in tracking data on the spacecraft’s trajectory or 
orbit. This program has also been utilized for the optimization of tracking schedules used during 
different phases of a space mission. 

Figures 4 and 5 show how the RMS errors in spacecraft position and velocity behave during 
the early phase of a transfer trajectory to the moon. The effects of different error sources are 
compared. On each graph, curve A shows the RMS errors in spacecraft position and velocity due 
to measurement noise in the tracking data. Curve B shows the effects of measurement noise and 
measurement bias, and Curve C shows the effects of measurement noise, measurement bias, and 
station location uncertainties. 




Figure 4— Position errors for the Apollo lunar 
transfer trajectory. 


Figure 5 — Velocity errors for Apollo lunar 
transfer trajectory. 


Because no a priori information is presumed about the state, the RMS errors in the state 
vector may be considered infinite. As soon as sufficient tracking data become available (a mini- 
mum of six independent measurements are required to define the components of the state vector), 
these RMS errors become finite and may be evaluated. 
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As increasing amounts of tracking data are utilized, the spacecrafts position and velocity er- 
rors are found to behave differently. Indeed, the RMS error in spacecraft position tends to reach 
an initial minimum rapidly and then increases monotonically, whereas the velocity continues to 
decrease monotonically. Several factors contribute to this difference in behavior: 

1. As the spacecraft’s distance from earth increases, the spacecraft - tracker geometry is 
weakened, resulting in an increase in the RMS errors in spacecraft position. 

2. The motion of a spacecraft along a lunar transfer trajectory is directed away from perigee 
towards apogee, resulting in a decrease in velocity (since maximum velocity is at perigee and 
minimum at apogee). The effect of decreasing velocity coupled with increasing distance from 

the earth causes a corresponding decrease in the RMS error in spacecraft velocity and a corres- 
ponding increase in the RMS error in spacecraft position. 

3. The RMS errors, when evaluated at discrete points along the spacecraft’s trajectory, are 
influenced by local anomalies such as the strong dependence on where along the trajectory these 
errors are evaluated, the number of tracking stations which track the spacecraft, etc. If the RMS 
errors in spacecraft position and velocity were evaluated at some fixed targeting point such as 
arrival at the moon’s region of dominance, then RMS errors in both position and velocity would 
tend to decrease as increasing amounts of tracking data become available. The needs of individual 
missions determine whether to evaluate the RMS errors at discrete points or at a fixed targeting 
point. For the Apollo mission, the instantaneous knowledge of the state (simulated by evaluating 
the RMS errors at discrete points) is required for updating the onboard computer. Likewise, for 
this mission, knowledge of the RMS errors evaluated at some targeting point (the lunar region of 
dominance) provides a basis for any corrective space maneuvers which may be required to mini- 
mize the errors at the targeting point. For the error analysis study described here, the RMS 
errors in spacecraft position and velocity are evaluated only at discrete points along the trajectory. 

Figures 4 and 5 show that the measurement bias errors are the principal contributors to the 
RMS errors in spacecraft position and velocity. The uncertainty in tracking station location be- 
comes less significant as the spacecraft’s distance from the earth increases. For the error 
analysis study described, the bias errors are not solved for but their effects are included; this 
results in a conservative estimate of the RMS errors in the state vector. Solving for some of 
the bias errors would provide a more optimistic estimate. However, to insure that a safety fac- 
tor is associated with the results presented in the error analysis, it is more prudent not to solve 
for any of the bias errors. 
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Appendix a 


Fundamental Matrices in Expanded Form as Used in This Report 


A listing of matrices used is presented here from Part I for completeness and convenience. 


SK (3X1) S P(3X1) 


SK (3X1) ( 3X1 ) 


(3X1) 


(3X1) 


(a-1) 


(a-2) 


sx 


(3X1) 


Sx l 

Sx 2 
Sx 

L (3xl ) 


(a- 3) 


SX 


(3X1 ) 


Si, 

Sx, 

Sx, 


(3X1 ) 


(a-4) 




(3x3) 


cose 0 sine 
0 10 
-sine 0 cos e 


_l(3*3) 


(a- 5) 


[R,(7t/ 2 - a)] 


(3x3) 


sin a cos a 0 
-cos a sin a 0 
0 0 1 


(3^3) 


(a-6) 


J(3*3) s [ R 2 <- s ) R 3 (w/2-a)J 


(3*3) 


(a-7) 
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ll ■ l 



(3*3) 


"1 0 0 “ 

0 - r cos e 0 

0 0 r 


-1(3*3) 


(a- 8) 


D° 


‘(1x3) 


3 [1. 


0] 


(1*3) 


(a-9) 


D? 


"(2x 3 ) 


0 -r cos e 
0 0 


cr 


r 


(2*3) 


(a- 10) 


and 


V 11 

V 12 

V 13~ 

V 12 

V 22 

V 23 

V 1 3 

" V 23 

V 33_ 


(a- 11) 


where 


V 


11 


= 0 


1 /• • . 

= - — (z t cos a - z 2 sin a) ) 


13 “ e = - ^-(z Jl sin £ sin a + z 2 sin £ cos a - z 3 cos £) 


( • r \ - S 0 c £ • 

£ tan £ - — ■ 1 = (z j sin a + z 2 cos a) 


-tan £ 


(zi cos a - z 2 sin a.) 


v 
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-1 


(Zj cos £ sin a + z 2 cos £ cos a + £ 


sin £ ) 


v 


12 


v 


21 ’ 
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and 


V 13 = " V 3 


Continuing the listing of matrices, 


[R t ( tt /2 - <P)] (3x3) - 


0 0 


0 s in cp cos (p 
0 -cos s in <p 


— \3 *3) 


(a-12) 


and 


where 


[R3U/2 + \>] (3x3) 


'-sin A. cos K 0' 
-cos K -sin K 0 
0 0 1 


C3><3) 


(a- 13) 


R (3«3) S [Rj <W/2 - 9> R s (V2+\)] 


(3X3) 


(a-14) 


^ (3x3) 


cos 6- sin (L 0~ 

G G 

-s in 0_ cos 0 

G G 


>(3x3) 


(a-15) 




(3-3) 


(a- 16) 


d t 


7-29211585 x 10~ s rad. /sec. 


&e r 


{«3<*c>} 


(3x3) 


-sin 0- 


cos 0 

G 


-cos (L -s in 0 

G G 


0 0j ( 


’(3x3) 
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Continuing the listing of matrices 


(a- 17) 


L(3*3) S [J T RR 3 (^)] 


(3*3) 


“(3-3) S [^ T R83^)^{J TR [^ R 3^o)]} ] 


(3x3) 


F, [l, 0 , ^ ] ( 1 x 3 ) 


‘(1X3) 


(2*3) 


"0 1 0” 
0 0 1 


(2*3) 


L °r "F F r L l 

r ( 1 * 3 ) L J 


(1*3) 


L° s Tf Li 

3 ( 2 * 3 ) L J ( j x 3 ) 


■ [ F ’"J 


(1*3) 


M? 


(2x3) 


[f.»] 


(2*3) 


M 3 (-V\ 3 , 3) 


"cos Q s -sin fi s 0 
sin Q g cos 0 

0 0 1 


(3*3) 


[Rj (- i s ) ] 


(3*3) 


10 0 
0 cos i g -sin i g 
0 s in i cos i 


(3*3) 


[R 3 (- w s )] (3x3) 


cos co -sin oj 0 

s s 

sin w cos w 0 

s s 

0 0 1 


—1(3-3) 


(a-18) 
(a- 19) 

(a-20) 

(a-21) 

(a-22) 

(a-23) 

(a-24) 

(a-25) 

(a-26) 

(a-27) 
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(a-28) 


■ ( 3.3) s [B 3 C-n.)«i<-i.)H3 (-«.>] 


(3x3) 


m -w 

L J (3^3) 


(3x3) 


-s in fi -cos n 0 ‘ 

s s 

cos Q. -sin H 0 

s s 

0 0 0 


Ri(-i s )R 3 (-s) 


'(3X3) 


r 3R ' 

l di s 

•N, 

L- a - 

(3X3) 


(3x3) 


r 

-o 

0 

0 

1 


0 

-sin i 

-cos i 

1 

) R 3 (- n s) 


s 

s 

R 3 (-"s) ( 

l 

0 

cos i 

5 

-s in i 

s ___ 

J 


(3x 3 ) 


BR 

B 


Kl M R 3 (- n s)Ri(- i s) 


(3x3) L J ( 3X3 > 


-sin a) 

s 

COS W s 

0 


-cos co 0' 
s 

-sin co 0 


— 1 S (3x3) 


; (3xi) 


a s ( COS E s - e S ) 

a s (1 - e s ) 1/2 sinE s 
0 


—1(3*1) 


and 


U 


(3*1) 


-n s a s sin E s 
(1 -e s cos E s ) 

n s a s (l~ e s) 1/2 cos E s 

(1 -e s cos E s ) 


0 


(3*1) 


where 


= (M ,/2 a -3/2) 


(a-29) 


(a-30) 


(a-31) 


(a-32) 


(a-33) 
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(a-34) 


Continuing the listing of matrices, 


[*q, * 1*.."'!*”."] 


(3X3) 


and 


"<.■»■ [“a, u !“i. "!>."] 


(3X3) 


where 


(3 x 3 ) 


m ll 

m 12 

m i 3 





m 2 1 

m 22 

"23 

"si 

" 3 2 

"33 


( 3 X 3 ) 


(cos E s - e s ) + 


3n s (t -r) sin E g 


2 (1 -e s cos E s ) 
sin 2 E r 


]• 


1 + 


(1-e cos E, 


) _ ’ 


a n sin E 
s s s 


21 


(l-e cos E s ) 


(1 - e 2) 1/2 sin E s 


3n s (t -r) (1 - e 2 ) 1/2 _cos E s 
2 (1 -e cos E s ) 


and 


a sin E s (cos E s - e g ) 


(1 _ e 2 ) 1/2 (1 -e s cos E s ) 

~ a s n s C 1 -^) 172 COS E s 

(l-e s cos E s ) 


m 31 = m 32 =m 33 = 0 • 


Continuing the listing of matrices, 


( 3 * 3 ) 


"ll 

"12 

" 13 

"21 

"22 

"23 

"31 

"32 

"33 
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( 3 * 3 ) 


(a-35) 


(a-36) 


(a-37) 



where 


(1 - e s COS E s) 


3n s (t-T) (cosE s -e s ) 


~ n s 3 s sin E s 


(1 - e s COS E s) 


(cos E s -e s ) 

-- cos E + 

2 s (1 - e s cos E s ) 


= n s a s (cos E s- e s > 

(1 -e s cos E s ) 3 


-"sCi-es)^ 2 3n s (t-r) sin E s 

1 = cos E^ 

2(l- e s c ° s E s ) (l-e s cos E s ) 2 


(1 -e cos E s V 


)S E s (cos E s -e s ) (1 -e^) 1 / 2 sin 2 E s 
(1 _ e 2 )l /2 " (l-e s cosE s ) 


S S > 

2 \ 1/2 


n s a s C 1 ^) sinE S 
(1 - e cos E s ) 3 


m 31 = m 32 = m 33 = 0 


Continuing the listing of matrices, 


[P(t)L 


< 3 * 3 > (6x 6) 


J 


*11 

*12 

*13 

*14 

*15 

*16 

*21 

*22 

*23 

*24 

*25 

*26 

*31 

*32 

^33 

*34 

^35 

*36 

*41 

*42 

^43 

*44 

^45 

*46 

*51 

*52 

^53 

*54 

^55 

*56 

*61 

*62 

^63 

*64 

^65 

*66 
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where 


0 n = cos ^ cos ^ » 

^12 = ^13 = ^14 = 0 , 

0 15 = - r s cos S sinX' t 

0, „ = - r s in S cos X. ' , 

^16 s ’ 

0 2 i - cos S sin X # , 

*22 = *23 = *24 = 0 ■ 

0 2S = r g cos 8 cos^' , 

0 O , = - r sin 8 sin X 1 . 

r 26 s 

^31 = sin S , 

*32 = *33 = *34 = *35 = 0 • 

*36 = r s COS S ’ 

*41 = 0 • 

0 42 = [siny s cos 8 cos V - cos y g sin a. sin V - cos y s cos sin 8 cos V] , 

0-, = v [cosy cos 8 cos X' + sin y sin a. sin V + siny cos a. sin S cos V] , 

s s si si 

0 44 = v s [- cos y s cos cl sin V + cosy s sin a. sin 8 cos V] , 

0. c = v [- siny cos S sin X' - cos y sin a. cos V + cosy cos a. sin 8 sin X' ] , 

0- ^ = v [- siny sin 8 cos V - cos y cos’ a. > cos S cos V] , 

r 46s 's ' s i J ’ 

*si =0 • 

0 C _ = [siny cos 8 sin A/ + cosy sin a. cos X 1 - cosy cos a. sin S sin X'] 

oz s SI SI * 

0 C , = v [cos y cos 8 sin X' - sin y sin a. cos X' + sin y cos a. sin S sin X # ] 

jo s s si si 

0 54 = v s [cosy s cos cl cos X' + cos y g sin a. sin 8 sin X'] , 

0_ c = v [sin y cos 8 cos X' - cos y sin a. sin X' - cos y cos a. sin £ cos X'] 

5 D S S SI S X 

0 CC = v [- siny sin 8 sin X' - cos y cos a. cos 8 sin X'l , 

* 6 , = 0 • 

*62 = [sin y s sin S + cos y s cos a. cos 8] , 
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0 63 = v g [cosy s sin 8 - sin y s cos a i cos S] # 
0,. = -v cosy sin a. cos S . 

r 64 s 's 1 » 

^ 65 -° > 

and 

0— = v [siny cos 5 - cosy cos a. sin 8] . 

OO S S S 1 
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Appendix b 


Coordinate Systems and Their Transformations 

This appendix lists the equations which relate the inertial coordinate system to the local co- 
ordinate system. All the transformations considered here, previously given in Part I, are renum- 
bered to be consistent with the order of presentation of the coordinate systems in Part II. 

1. The inertial coordinate system is a right handed Cartesian coordinate system defined as: 
(Figure 1, Part 1) 


r*; 


(3*1) 


(3*1) 


(b-1) 


where 

Origin is: Earth f s center of mass 
x x -axis: Directed towards vernal equinox. 
x 2 - axis: Normal to both x 1 and x 3 -axes. 
x 3 -axis: Directed along earth's axis of rotation. 

2. The local coordinate system is a right handed Cartesian coordinate system defined as: 
(Figure 1 and 3, Part I) 


(3*1) 


(3*1) 


(b-2) 


where 

Origin is: At the observer 

Zj-axis: Directed towards local east. 

z 2 -axis: Directed towards local north. 

z 3 -axis: Directed along normal to local horizon plane. 
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3. The coordinates of the observer. (Figure 4, Part I) 



" S l _ 


"(N +h) 

cos tp cos A. 

S (3xl)~ 

S 2 

= 

(N +h) 

cos c p sin A. 


_ S 3_ 


[N(l - 

e 2 ) + h] sin cp 


Origin: The center of the reference ellipsoid 


(b-3) 


and 


N = a 0 /(l - e 2 sin 2 cp) 1/2 


4. Transformation equations relating the local coordinate system with the inertial coordinate 
system (Figure 1, Part I) 

^(3x1) " ^(3x3) ®3^(3x3) ^(3x1) " ®(3x3) ^(3x1) * (k-4) 

where 


R (3*3) = [R 1 (7r/2 - <P) ®3 < w/2 + ^ ( 3 « 3 ) 

(see Appendix a for the expanded form of matrices R (3x3) and R 3 (0 G )). 



and 


0 < 0 G < 2tt , 


or 


= Greenwich sidereal time 

G 

<9 = Greenwich sidereal time at 0 h U.T. 

G o 

= 15?041 06864/h r 



= 7.29211585 x 1CT 5 radians/sec. 


t : Universal time of observation. 

5. The transformation from the local coordinate system into the coordinates of the spacecraft 
in its orbit is obtained by combining equation (c-20) with (b-4): 

Z ( 3 xi)= U (3xi) - R (3*3) S (3*i) . ( b “ 5 ) 

(3x3) 
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where 


R (3*3)^R3<- n >Rl<- « V-"»(3K3) 

and 


U, 


(3x1) 


a s (cos E s - e s ) 
a (1 - e 2 ) 1/2 sin E 

s x s' s 


0 





1 
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Appendix c 

Transformation of Position and Velocity Vectors 
into Keplerian Orbital Elements 


This appendix gives the transformation of the position and velocity vectors from the inertial 
coordinate system into the Keplerian orbital elements. Before doing so the transformation of the 
following spacecraft orbital injection parameters into the inertial coordinate system will be given. 

For the purposes of this paper, injection orbital parameters are defined as follows: 

P s — Radial distance of spacecraft from center of earth, 
v s — Magnitude of velocity vector. 
y s — Flight path angle. 
a . — Injection azimuth. 

A.' — Right ascension of spacecraft at time of injection into orbit. 

8 — Declination of spacecraft at time of injection into orbit, 
r — Epoch time (time spacecraft is inserted into orbit). 

The transformation of the injection orbital parameters into the position and velocity vectorial 
components in the inertial coordinate system are: 



" x l~ 


P s cos S cos k‘ 

p* = 

X 2 

= 

p cos S sin k * 

1 s 


_ X 3 _ 

(3x1) 

p sin S 


(c-l) 


X 


1 


(3x1) 


v s {siny s cos S cos k 1 - cosy s sin a. sin k' 
v $ {siny s cos S sin k* + cos y s sin a. cos A.' 
v s {siny s sin S + co s y g cos a. cos 8 } 


- cosy s cos a. sin 8 cos A/} 

- cosy cos a. sin 8 sin A.'} 

' s 1 


(c-2) 


1 (3x1, 


From the position and velocity vectors given by (c-l) and (c-2), the Keplerian orbital elements 
are derived. The geometry of the conditions at the time of insertion into orbit is shown in 
Figure c-l. 
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For the determination of the orbit orien- 
^ ^ directed along earth's tation with respect to the inertial coordinate 

ROTATIONAL AXIS 

system, the following vectors are required: 

SATELLITE ORBIT 



(VERNAL EQUINOX) 


Figure C-l — Geometry of satellite orbit in the 
inertial reference frame. 


v° 


n° = 


s 

1 - * 1 

(c-3) 



' IvJ 

(c-4) 

? s x v s 

(c-5) 

x V.l 


. sin e Q (n ° x 

(c-6) 

cos 0 O (n ° x p°) 

(c-7) 

n° x k° 


rf° x k°| 

(c-8) 


From the energy integral the semi- major axis is determined: 


fj_ _ {ivj^ v 1 
\Mi m / 


(c-9) 


The eccentricity is given by: 


e . = M . {|p,xv.l>^ 1/a 


(C-10) 


The true anomaly is determined from the polar form of the equation of an ellipse. 


cos 6 q = . 


a (1 - e 2 ) 

s ' « 7 

IpJ 


-i 


sin e 0 = 


.(1 - el) 


2AI/ 2 


\(P°s ’ 7 s)f 


(c-11) 
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where & 0 - O s (r). From the true anomaly the transformation into the eccentric anomaly at epoch 
is given by: 


E° = 2 tan' 


■i ( - - c ° s g °y /2 i 

Lx 1 + e J V + cos 0 o) J 


The mean anomaly at epoch is then derived from Kepler 1 s equation: 


M o = < E :- e s sinI 0 


To evaluate the eccentric anomaly at times other than at epoch, let 


(c-12) 


(c- 13) 


M s = M 0 + 0* 1/2 a" 3/2 ) (tj - r) 


(c-14) 


where ( tj - r)is to be expressed in seconds of time. 

M 0 is in radians, 
r is epoch time. 

Then: 


E<°> = M 

S S 

= M 

5 S 

E< 2 > = M 

s s 

E<"> = M + e sin E (n_1 > 

s s s s 


+ e sin 

s s 

+ e sin E^ 1 ^ 

s s 


(c-15) 


That value of E s is then chosen for which the following convergence criterion is met, 


{E< n > - (M s + e s sinE< n_1 ))} < l(T 8 


(e- 16) 


The orbital plane orientation remains invariant whenever a Kepler ian -type orbit is pre- 
sumed. Using the fundamental vectors given by (c-3) to (c-8), the following orientation param- 
eters are obtained. 
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Inclination of orbit to equatorial plane: 


i g = cos” 1 (n° • k°) 


0 < i < 77/2 . 

s 


(c-17) 


Argument of perigee: 


sin a) 

s 


cos 6 0 - sin 6 q (V • p*) 


cos w = ( / t° m p°) 
0 < oj < 2tt . 

s 


(c-18) 


The right ascension of the ascending node: 


5 i n q _ a°-p?) - 

s (p;-j°) 


cos n s = (V- r°) 

0 < n < 277 


(c— 19) 


Once the position and velocity vectors at epoch have been transformed into the orbital ele- 
ments, it is possible to generate an ephemeris of the satellite’s motion in the Keplerian sense 
using the following transformation equations. 

The position vector at any time t is given by: 

^(3x1) = ^3 ^3^ W s^(3x3) ^(3x1) • (C”20) 


where 


P s (t) ~ 


Xj ( t) n 

x 2 (t) 


x 3 (t) 


(3x1) 
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(c— 21) 


U 


(3x1) 


a # (cos E s - e s ) 

a s ^ - e s) 1/2 sin E s 

0 


E s = M s + e s sin E„; M s = M 0 + O 1 ' 2 a' 3/2 ) (t - r) , 


s “ ‘”0 


and the velocity is given by 


^(3*1) “ ®3 ^ O ^3 W s^(3x3) ^(3x1) 


where 


v.(t) 


(3x1) 


Xj(t) 

x 2 (t) 


U 


(3xl) ” 


- n a sin E 

s s s 

(1 - e cos E ) 

v s s' 

n a (1 - e 2 ) 1/2 cos E 

s s ' s s 

(1 - e s cos E s ) 


0 


_l (3x1) 


n = u} /2 a~ 3/2 ; n = Mean motion , 

s “ s s ’ 


and 


ul, ., s = 3.986032 x 10 s km 3 /sec 2 . 

“(earth) 

By assuming no perturbations, a Keplerian orbit adequately describes the spacecraft’s motion 
around the earth. Utilizing Equations (c-20) and (c-21) the position and velocity can be calculated 
for any instant of time. 


NASA- Langley, 1966 
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